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ABSTRACT 


The  three-dimensional  flow  in  the  region  where  a  turbulent  boundary  layer  encounters  an 
appendage  is  characterized  by  the  formation  of  a  junction  root  vortex.  The  formation  of  this 
vortex  can  be  attributed  to  the  adverse  pressure  gradient  upstream  of  the  appendage.  One 
major  area  of  concern  is  the  temporal  behavior  of  the  vortex.  Fluctuations  in  the  size  and/or 
location  of  the  vortex  can  cause  degradation  of  downstream  components  or  excitation  of 
local  structures.  Recent  studies  in  wind  tunnels  have  identified  a  bi-modal  characteristic  in 
the  structure  of  the  vortex. 

An  experimental  investigation  of  the  time-dependent  behavior  of  the  junction  root  vortex  was 
conducted  using  MIT’s  variable  pressure  water  tunnel  and  laser  doppler  anemometry  system. 
Data  was  collected  at  various  points  forward  of  a  modified  NACA  0020  section  appendage 
mounted  on  a  splitter  plate  and  analyzed  using  off-line  DFT  analysis.  The  data  was  gathered 
at  tunnel  speeds  of  10  and  20  ft/sec  using  an  analog-to-digital  acquisition  system.  Velocity 
probability  distributions  using  a  fully  digital  acquisition  system  were  also  plotted  to  identify 
the  bi-modal  structure. 

The  spectra  generated  from  the  DFT  analysis  were  inconclusive.  No  distinct  frequencies  or 
pattern  of  frequencies  attributable  to  the  flow  could  be  identified.  Fluctuations  within  the 
spectra  made  their  accuracy  questionable.  The  existence  of  the  bi-modal  structure  was  con¬ 
firmed  by  the  velocity  probability  distributions.  -r  . 
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1  INTRODUCTION 


1.1  BACKGROUND 

The  two-dimensional  turbulent  boundary  layer  along  a  body  has  been  studied  exten¬ 
sively  and  is  well  understood.  When  this  boundary  layer  encounters  an  appendage,  how¬ 
ever,  the  flow  in  the  region  of  the  junction  becomes  extremely  complex.  This 
three-dimensional  flow  is  characterized  primarily  by  the  formation  of  a  "horseshoe"  or 
junction  root  vortex  at  the  base  of  the  appendage.  There  are  many  practical  examples  of 
this  flow  including  the  keel  of  a  sailboat,  a  submarine’s  sail  and  control  surfaces,  an  air¬ 
plane’s  wings  and  stabilizers,  and  even  a  telephone  pole  in  a  windstorm. 

The  formation  of  this  vortex  secondary  flow  can  be  attributed  to  the  adverse  pres¬ 
sure  gradient  upstream  of  the  appendage.  The  vorticity  in  the  incoming  turbulent  bound¬ 
ary  layer1  is  concentrated  or  coalesced  by  the  gradient  forming  the  vortex.  The  flow 
causes  the  vortex  to  wrap  around  the  appendage  turning  the  cross-flow  vorticity  of  the 
incoming  boundary  layer  into  streamwise  vorticity.  The  vortex  produces  a  three- 
dimensional  separation  region  in  the  vicinity  of  the  appendage. 

Dickinson  [5]  divided  the  flow  around  an  appendage  into  three  conceptual  regions. 
The  inflow  region  is  the  region  upstream  of  the  appendage  where  the  presence  of  the 
appendage  is  not  felt.  The  second  region  is  the  appendage  region.  This  is  the  region 
which  includes  the  adverse  pressure  gradients  and  the  generation  of  the  secondary  flows 
associated  with  the  appendage.  This  region  is  further  divided  into  the  leading  edge  subre¬ 
gion ,  the  alongside  subregion ,  and  the  trailing  edge  subregion.  The  final  region  is  the 
wake. 


1  A  junction  vortex  can  also  form  when  the  incoming  boundary  is  laminar. 
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A  numbe  of  experimental  studies  have  been  conducted  or  are  currently  in  progress 
to  characterize  the  junction  vortex.  Most  of  the  early  work  concentrated  on  investigating 
the  flow  for  various  appendage  geometries.  Love  [9]  conducted  some  preliminary 
research  which  showed  that  the  curvature  of  the  leading  edge  affected  the  strength  of  the 
generated  vortex.  Since  then,  Belik  [2]  and  Baker  [1]  have  investigated  the  flow  using 
cylindrical  appendages;  while  Shabaka  and  Bradshaw  [21],  Oguz  [13],  Kubendran  [8], 
and  McMahon  et  al.  [10], [1 1]  conducted  studies  using  infinite-chord  appendages.  To 
study  the  flow  around  more  realistic  hydrodynamic  shapes,  Dickinson  [4], [5]  of  the 
David  Taylor  Research  Center  (DTRC)  and  Pierce  et  al.  [14], [15], [16]  of  Virginia  Poly¬ 
technic  Institute  (VPI)  have  used  streamlined  foils  in  their  experimental  work.  A  more 
complete  literature  review  on  the  topic  of  the  junction  vortex  can  be  found  in  Sung  [22], 
This  document  not  only  references  experimental  studies,  but  also  publications  on  analytic 
and  numerical  work. 

The  majority  of  experimental  work  on  the  junction  vortex  flow  has  been  for  the  pur¬ 
pose  of  obtaining  data  to  support  the  development  of  analytic  or  numerical  prediction 
schemes.  The  data  has  been  collected  using  various  experimental  techniques  including 
flow  visualization,  hot -wire/film  anemometry,  laser  doppler  anemometry,  and  pressure 
taps  with  most  of  the  investigations  having  been  conducted  in  low-turbulence  wind  tun¬ 
nels.  Recently,  a  great  deal  of  attention  has  been  focused  on  reducing  the  strength  of  the 
generated  root  vortex.  In  marine  applications,  the  presence  of  this  vortex  can  degrade  the 
performance  of  downstream  control  surfaces  and  propellers.  The  use  of  fillets  and  other 
leading  edge  modifications  is  currently  under  study  by  various  institutions.  Pierce,  Fran- 
gistas,  and  Nelson  [17]  have  published  some  preliminary  results  from  the  studies  being 
conducted  at  VPI. 
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Another  area  of  concern  is  the  temporal  behavior  of  the  junction  vortex.  Fluctua¬ 
tions  in  the  size  and/or  location  of  the  vortex  can  cause  additional  degradation  of  down¬ 
stream  components,  or  excitation  of  local  structures  resulting  in  vibrational  and  acoustical 
problems.  Rood  et  al.  [18], [19]  conducted  some  of  the  first  studies  on  the  time-dependent 
nature  of  the  junction  vortex  flow.  He  reported  relatively  low  frequency  and  narrow' 
bandwidth  turbulent  fluctuations  in  the  three-dimensional  boundary  layer  downstream  of 
a  streamlined  protuberance  with  a  1.5:1  elliptical  nose,  similar  to  the  one  used  in  this 
study.  In  a  meeting  on  appendage  design  at  DTRC  in  1988,  Dr.  R.L.  Simpson  [3]  of  VPI 
reported  detecting  a  bi-modal  characteristic  in  the  structure  of  the  vortex.  Histograms  of 
streamwise  velocity  fluctuations  measured  close  to  the  line  of  shear2  exhibit  two  distinct 
peaks.  Pierce  et  al.  [17]  have  shown  that  the  spectra  within  the  vortex  flow  is  similar  to 
that  in  the  upstream  boundary  layer,  suggesting  that  the  fluctuations  in  the  large  scale  vor¬ 
tex  structure  originate  in  the  boundary  layer. 

1.2  RESEARCH  GOALS 

This  experimental  investigation  had  two  purposes.  The  primary  goal  was  to  gather 
raw  data  on  the  junction  vortex  flow  and  process  it  to  obtain  information  on  the  time 
dependent  nature  of  the  flow  through  spectral  analysis.  Unlike  previous  studies  which 
were  conducted  using  wind  tunnels,  this  investigation  was  conducted  in  a  circulating 
water  tunnel.  An  attempt  to  verify  the  bi-modal  structure  identified  by  Simpson  and 
quantify  its  frequency  of  fluctuation  was  a  principal  element  of  this  objective. 


2  The  line  of  shear  delineates  a  region  of  high  shear  near  the  foil  from  an  area  of  lower  shear  outside.  The  line 
of  shear  falls  inside  the  separation  line  forward  of  the  appendage. 
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The  other  purpose  of  the  investigation  was  to  determine  the  feasibility  of  using 
velocity  data  from  the  laser  doppler  system  in  MIT's  variable  pressure  water  tunnel  facil¬ 
ity  to  conduct  spectral  analysis.  A  key  element  of  this  goal  was  to  incorporate  a  new 
all-digital  acquisition  and  flow  analysis  system  into  the  laboratory  facility. 
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2  EXPERIMENTAL  APPARATUS 

2.1  VARIABLE  PRESSURE  WATER  TUNNEL 

All  testing  and  data  acquisition  was  conducted  using  MIT's  variable  pressure  water 
tunnel  [7].  This  facility  is  designed  to  test  and  evaluate  the  dynamic  performance  of  pro¬ 
pellers  and  two  and  three  dimensional  hydrofoils.  The  test  section  is  a  54  inch  long,  20 
inch  square  closed-jet  section.  The  maximum  test  section  speed  is  33  ft/sec. 

The  area  directly  upstream  of  the  test  section  contains  the  upstream  turning  vanes,  a 
circular  to  square  transition  section,  a  settling  section,  and  a  contraction.  The  transition 
section  has  a  contraction  ratio  of  1.46:1,  while  the  ratio  for  the  contraction  section  is 
4.92:1.  The  overall  contraction  ratio  for  the  facility  is  7.17:1.  The  settling  section  con¬ 
tains  a  honeycomb  consisting  of  one  inch  diameter  acrylic  tubes,  each  one  being  one  foot 
long,  and  three  wire  screens.  The  settling  section  is  designed  to  reduce  the  turbulence  in 
the  test  section.  Measurements  of  the  free  stream  turbulence  using  the  laser  doppler  velo- 
cimeter  show  a  turbulence  intensity  of  less  than  0  5  percent,  accounting  for  phase  noise3 
in  the  velocity  measurements. 

The  test  section  is  immediately  followed  by  a  diffuser,  downstream  turning  vanes, 
and  a  square  to  circuL-r  transition  section.  The  closed  circuit  flow  in  the  water  tunnel  is 
driven  by  a  30.25  inch  diameter  four-bladed  impeller. 

2.2  BODY/APPENDAGE  JUNCTION  MODEL 

The  model  used  in  the  study  of  the  junction  vortex  consisted  of  a  modified  NACA 
0020  foil  section  mounted  on  a  flat  splitter  plate.  The  foil  and  plate  were  both 


3  Phase  noise  is  related  to  the  bandwidth  of  the  doppler  signal  spectrum  and  can  be  related  to  an  equivalent 
turbulence  intensity. 
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constructed  using  anodized  aluminum.  The  foil  section  was  produced  on  a  numerically 
controlled  milling  machine. 

Figures  1  and  24  show  the  model  layout,  test  section  mounting  arrangement,  and 
experimental  coordinate  system. 


Figure  1. 

Test  Section  Layout 

(From  View) 


4  All  dimensions  in  Figures  1  and  2  are  in  inches. 
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Figure  2. 

Test  Section  Layout 

(Side  View) 

The  shape  of  the  foil  section  was  chosen  to  conform  to  several  previous  investiga¬ 
tions  (see  Section  1.1).  The  standard  NACA  0020  section  was  modified  by  replacing  the 
nose  with  a  1.5:1  elliptical  arc.  This  bluff  leading  edge  produces  a  stronger  adverse 
pressure  gradient  which  promotes  the  formation  of  the  junction  vortex.  Equations  2. 1  and 
2.2  were  used  to  generate  the  foil  section  offsets. 

Nose  Section: 

v2  (x-1.8)2  ,  t2-1) 

—.  + - ; — =  1 

1.22  1.82 

Tail  Section: 

y/c  =  0.2969  ■  Vx/c  -0.126  •  x/c  -  0.3516  •  (x/c)2  (2-2) 

+0.2843  •  (x/c  )3  -  0. 1 0 1 5  •  (x  'c  )4 
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The  two  shapes  were  joined  at  their  respective  maximum  widths  (2.4  inches) 
resulting  in  a  chord  length  of  10.3  inches.  The  span  of  the  foil  was  set  at  15  inches  to 
allowing  a  2.4  inch  gap  between  the  foil  and  top  of  the  test  section. 

6 

5  - 

4  - 

3  - 
2  - 


Z 


0  2  4  S  8  10  12 

CHORD  LEIMQTH  (IN.) 


Figure  3. 

Modified  NACA  0020  Foil  Section 

Trip  wires  of  circular  cross-section  were  positioned  just  aft  of  the  leading  edge  of 
the  splitter  plate  and  0.4  inches  upstream  of  the  maximum  thickness  of  the  foil.  The  trip 
wire  on  the  splitter  plate  invoked  the  development  of  a  turbulent  boundary  layer  along  the 
plate  forward  of  the  foil.  As  recommended  by  Devenport  and  Simpson  [3],  the  boundary 
layer  on  the  foil  was  tripped  to  avoid  any  unsteadiness  of  the  flow  which  might  have 
resulted  from  natural  transition.  A  trip  wire  diameter  of  0.047  (3/64)  inches  was  used  in 
both  cases  making  them  fully  effective  down  to  a  tunnel  speed  of  2.26  ft/sec  as  stated  by 
Schlichting  [20]. 
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2.3  LASER  DOPPLER  VELOCIMETER  SYSTEM 


The  velocity  data  used  for  this  study  was  gathered  using  the  horizontal  component 
of  a  three-component  laser  doppler  velocimeter.  A  three  watt  Argon-Ion  laser  serves  as 
the  coherent  light  source  for  the  measurement  system.  The  beam  of  light  from  the  laser  is 
split  into  three  components,  green  (X=5.14  nm),  blue  (X=488.0  nm),  and  violet  (^=476.5 
nm),  by  an  optical  color  separator.  The  blue  component  is  used  for  horizontal  measure¬ 
ments. 

Each  component  of  the  laser  is  further  split  into  two  beams  in  a  dual  beam  system. 
By  use  of  optical  lenses  and  mirrors,  the  two  beams  intersect  at  a  fixed  angle  inside  the 
tunnel  test  section  forming  an  ellipsoidal  measurement  volume.  The  two  intersecting 
beams  of  light  produce  a  fringe  pattern  within  this  measurement  volume  perpendicular  to 
the  plane  formed  by  the  two  beams.  A  particle  in  the  fluid  moving  through  this  fringe 
pattern  will  scatter  light.  The  intensity  of  the  scattered  light  will  be  modulated  at  a  fre¬ 
quency  determined  by  the  fringe  spacing  and  velocity  of  the  particle  perpendicular  to  the 
fringe  pattern.  This  frequency  is  called  the  doppler  frequency  and  is  given  by  Equation 
2.3. 


2  m  sin  k 

X 


(2.3) 


fD  =  doppler  frequency 

u  =  component  of  velocity  perpendicular  to  the  fringe  pattern 
k  =  one-half  of  the  intersection  angle  of  the  two  beams 
X  =  wavelength  of  light 


The  scattered  light  from  particles  crossing  the  measurement  volume  is  collected  by 
a  set  of  receiving  optics  and  detected  with  a  photomultiplier  tube  which  converts  the 
modulated  light  signal  to  a  voltage  signal  w  ith  the  same  doppler  frequency.  The  signal  is 
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then  processed  by  a  data  acquisition  system  to  obtain  the  measured  component  of  the  par¬ 
ticle  velocity.  The  receiving  optics  can  be  placed  at  any  location  where  the  measuring 
volume  can  be  effectively  viewed.  In  this  study,  the  same  optics  were  used  to  focus  the 
outgoing  light  and  collect  the  scattered  light;  this  is  called  an  on-axis  backscatter  system. 

A  stationary  fringe  pattern  cannot  differentiate  between  equal  positive  and  negative 
flow  velocities,  nor  can  it  convey  a  zero  flow  condition.  To  overcome  this  shortcoming,  a 
frequency  shift  is  applied  to  one  of  the  two  crossing  beams.  This  causes  the  fringe  pattern 
to  slowly  move  upstream  or  downstream  depending  on  the  magnitude  and  direction  of  the 
shift.  The  moving  fringe  pattern  will  produce  a  doppler  frequency  even  in  a  zero  flow 
condition  and  will  provide  different  frequencies  for  equal  positive  and  negative  flow 
velocities. 

The  laser  doppler  velocimeter  used  in  this  study  was  mounted  on  a  three-axes  tra¬ 
verse  table.  The  motion  of  the  table  was  controlled  remotely  from  a  personal  computer 
using  feedback  control.  The  positional  accuracy  of  the  traverse  table,  as  stated  in  the 
operations  manual,  is  2.5  pm. 

Additional  information  on  laser  doppler  anemometry  theory  and  techniques  can  be 
found  in  Durst,  Melling,  and  Whitelaw  [6]. 

Table  I. 

Laser  System  Characteristics 


Wavelength  488.0  nm 


Diameter  of  63  pm 

Measurement  Volume 


Length  of  736  pm 

Measurement  Volume 


Number  of  Fringes  22 


Fringe  Spacing  2.85  pm 
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3  DATA  ACQUISITION 


j 

3.1  ANALOG-TO-DIGITAL  ACQUISITION 


The  data  for  the  spectra  in  Appendix  B  was  collected  using  an  analog-to-digital 


acquisition  system.  Figure  4  depicts  this  acquisition  system. 


Figure  4. 

Analog  Data  Acquisition  System 

As  shown,  the  modulated  voltage  signal  from  the  photomultiplier  tube  was  fed  to  an 
analog  tracker.  The  tracker  is  a  signal  processor  specifically  designed  to  convert  the 
doppler  frequency  to  an  analog  voltage,  and  extract  the  signal  from  the  background  noise. 
Since  the  signal  from  the  photomultiplier  tube  is  discontinuous  and  varies  in  amplitude. 
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the  tracker  serves  as  a  combination  of  tracking  filter,  frequency  to  voltage  converter,  and 
sample  and  hold  circuit.  The  validity  of  the  signal  is  confirmed  before  making  it  avail¬ 
able  for  output  as  a  voltage  directly  proportional  to  flow  velocity. 

The  output  voltage  from  the  tracker  was  sent  through  a  signal  conditioner.  The  only 
purpose  of  this  device  was  to  apply  a  1kHz  low  pass  filter  to  the  voltage  signal,  which 
was  necessary  to  prevent  aliasing  during  the  off-line  spectral  analysis  (see  Section  4.1). 
The  filtered  voltage  was  then  converted  to  digital  information  by  using  one  channel  of  a 
16-channel  analog-to-digital  converter  board  installed  in  the  laboratory’s  personal  com¬ 
puter. 

The  operation  of  the  converter  board  was  controlled  via  software  directly  from  the 
computer  keyboard,  permitting  selection  of  sampling  frequency  and  number  of  data 
points.  The  digitized  voltage  versus  time  data  was  convened  to  an  output  file  and  placed 
on  a  floppy  disk  for  later  analysis.  For  each  collection  of  samples,  a  graphic  display  of 
voltage  versus  time  was  available  on  the  computer  screen.  Figure  5  is  an  example  of  a 
typical  display.  Note  the  large  amount  of  scatter  in  the  tracker  output. 
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Figure  5. 

Typical  Digitized  Data  Display 

The  digital-to-analog  converter  and  data  analysis  software  used  to  generate  the 
spectra  in  Appendix  B  was  designed  for  operation  on  continuous  signals.  The  tracker's 
output  signal,  however,  appears  as  "points"  or  "lines"  due  to  the  data  validation  circuity. 
This  simply  means  that  the  last  reading  is  held  until  a  new  reading  is  validated.  The  out¬ 
put  then  is  basically  a  series  of  point  measurements  rather  than  a  continuous  output.  To 
meet  the  continuity  requirements  for  the  analysis,  it  was  necessary,  therefore,  to  insure 
that  the  data  rate  from  the  tracker  sufficiently  exceeded  die  data  sampling  frequency.  In 
all  cases,  the  tracker  data  rate  exceeded  the  2000  Hz  sampling  frequency  by  a  factor  of 
five. 
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3.2  FULL  DIGITAL  ACQUISITION 


The  data  for  the  velocity  probability  distributions  in  Appendix  C  was  collected 
using  a  fully  digital  acquisition  system.  Figure  6  depicts  this  acquisition  system. 


Figure  6. 

Digital  Data  Acquisition  System 

The  primary  difference  between  the  data  acquisition  system  in  Section  3.1  and  that 
in  Figure  6  is  incorporation  of  a  digital  signal  processor  or  Intelligent  Flow  Analyzer 
(IFA)5.  One  exceptional  feature  of  the  IFA  is  that  it  makes  only  "good"  measurements  or 
it  makes  no  measurements  at  all.  The  IFA  determines  the  doppler  frequency  by  measur¬ 
ing  the  duration  of  eight  contiguous  cycles  of  the  doppler  signal.  This  elapsed  time  is 


5  The  Intelligent  Flow  Analyzer  (IFA)  is  a  product  of  Thermal  Systems.  Inc.  (TSI). 
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then  processed  to  obtain  the  frequency.  To  extract  the  true  doppler  signal  from 
background  noise,  the  IFA  uses  an  autocorrelation  of  zero-crossings  technique.  The 
result  is  highly  accurate  data. 

The  IFA  was  not  used  to  gather  the  data  for  the  spectral  analysis  for  two  reasons: 

(1)  It  was  not  available  during  the  initial  measurements.  (2)  The  processor  data  rate  was 
insufficient  to  meet  the  2000  Hz  sampling  requirement.  In  many  cases,  the  maximum 
obtainable  data  rate  was  lower  by  a  factor  of  ten  than  that  required  by  the  spectral  analy¬ 
sis  (see  Section  4.1).  Non-optimal  tunnel  seeding6  and  the  processors  stringent  validation 
algorithm  were  the  probable  causes  of  this  unfortunate  limitation. 

The  IFA  was  directly  controlled  from  the  laboratory  personal  computer.  The  real¬ 
time  display  software  installed  on  the  computer  not  only  controlled  the  IFA  functions,  but 
also  produced  the  velocity  probability  distributions  found  in  Appendix  C. 


6  Seeding  is  the  process  of  adding  particles  to  the  flow  to  improve  the  fluids  ability  to  scatter  light. 
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4  DATA  ANALYSIS 


4.1  SPECTRAL  ANALYSIS  FUNDAMENTALS 


The  purpose  of  spectral  analysis  is  to  determine  the  frequency  composition  of  a 
continuous  random  function  of  time  x(t).  Since  most  experimental  measurements  of  ran¬ 
dom  functions  are  now  conducted  digitally,  the  spectral  analysis  is  performed  on  a  dis¬ 
crete  series  of  numbers,  x/t),  obtained  by  sampling  the  continuous  function  x(r).  The 
introduction  of  the  fast  Fourier  transform  (FFT)  has  made  the  calculation  of  spectral 
estimates  directly  from  the  discrete  time  series  extremely  quick  and  very  accurate. 

Any  non-periodic  function  of  time  can  be  represented  by  a  Fourier  integral  of  the 

form 


x(t) 


•i 


i 


A  (w )  cos (wt  )dw  +  B(w)  sin(wr  )dw 


(4.1) 


w  =  angular  frequency  (rad/sec) 


where 


A(w) 


— r 

2rc  J-o 


x(t)cos(wt)dt 


1 

£(w)  =  —  x(t)sm(wt)dt 
2k 


(4.2) 


For  engineering  purposes,  Equation  4.1  is  true  in  classical  Fourier  analysis  if 


j  x(r)  |  dt  <  °° 


(4.3) 


Defining  X/w)  as 


X(w)  =  A(w)-iB{w) 


(4.4) 
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Equations  4.2  can  be  combined  to  give 


l  r~  ,  (4.5) 

X(w)  =  -J  x{t) e-’dt 

Equation  4.5  is  the  formal  definition  of  X(w)  and  is  called  the  Fourier  transform  of  x(t). 

The  condition  given  by  Equation  4.3  cannot  be  satisfied  by  an  unending  stationary 
process.  Therefore  the  classical  Fourier  theory  cannot  be  applied  to  most  experimental 
data.  This  difficulty  is  overcome  by  analyzing  the  autocorrelation  function  RJ z).  The 
autocorrelation  of  a  function  x(t)  is  defined  as  the  average  of  the  product  x(t)x(t+x),  which 
for  a  stationary  process  is  only  a  function  of  the  time  difference  t.  The  autocorrelation 
function  can  be  adjusted  to  meet  the  condition  of  Equation  4.3  allowing  it  to  be  Fourier 
transformed.  The  transform  of  RJ t)  is  called  the  spectral  density  or  the  spectrum  of  x(t). 

Most  experimental  data,  as  stated  above,  is  a  discrete  series  of  samples  xj t),  r  = 
0,1,2,. where  t  =  rA  and  A  =  T/N.  T  is  defined  as  the  total  sampling  time  and  N  is 
the  total  number  of  samples.  In  this  situa:'  jn,  the  integral  in  Equation  4.5  can  be  replaced 
by  the  summation 

w  1  —  l  — i2wt  r/N  iaiO  /\7  i\  ^  ^ 

Xk  — —2Lr _0  xre  k  —  0, 1,2,  1) 

Equation  4.6  is  the  formal  definition  of  the  discrete  Fourier  transform  (DFT),  where  k  is 
the  harmonic  number  that  is  related  to  a  discrete  angular  frequency  w  =  2nkiJ.  The  fast 
Fourier  transform  is  a  computer  algorithm  for  calculating  DFT's. 


7  A  process  is  said  to  be  stationary  if  the  probability  distribution  of  a  collection  of  samples  is  independent  of 
absolute  time. 
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The  advantage  of  the  DFT  is  that  it  can  be  used  to  generate  the  spectrum  of  x(t) 
directly  from  the  discrete  time  senes  xr(t)  without  first  generating  the  autocorrelation 
function.  The  theory  supporting  this  can  be  found  in  Newland  [12].  The  DFT  produces 
spectral  components  at  the  frequencies  defined  by  the  harmonic  numbers  k  = 

0,1. 2,..., (AM).  Values  of  the  spectrum  for  k  >  (N- 1)  can  be  shown  to  simply  repeat  them¬ 
selves.  Additionally,  it  can  be  shown  that  the  spectrum  is  symmetrical  about  the  zero 
frequency  position.  The  unique  and  valid  portion  of  the  spectrum  exists  only  between 
w  =  0  and  w  =  7t/A.  The  upper  limit  on  the  frequency  is  called  the  Nyquist  frequency8  vv0 
and  is  the  maximum  frequency  that  can  be  detected  from  data  sampled  at  a  time  spacing 
A.  If  frequencies  greater  than  w0  exist  the  original  signal,  distortion  of  the  spectrum 
called  aliasing  can  occur.  The  original  signal  must  be  filtered  prior  to  processing  to 
remove  these  frequencies  to  avoid  this  distortion. 

The  accuracy  of  a  spectral  measurement  depends  on  the  effective  bandwidth  of  ue 
measurement  and  the  total  sampling  time  T.  The  ratio  of  the  standard  deviation  a  to  the 
mean  value  m  of  the  spectral  density  is  given  by  Equation  4.7. 

£  =  _1__  (4.7) 

m  ~yjBJ 

For  most  practical  methods  of  calculating  the  spectrum  using  DFT’s.  the  effective 
bandwidth  B,  is  approximately  1/T  Hz.  This  means  that  aim  =  1 .  or  the  magnitude  of  the 
standard  deviation  of  the  spectrum  is  equal  to  the  mean.  To  improve  this  accuracy, 
several  adjacent  spectral  estimates  can  be  averaged  increasing  the  effective  bandwidth. 
Equation  4.7  becomes 


8  The  Ny  quist  frequency,  when  converted  to  Hz,  is /„  =  1/2A. 
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(2«+l)  =  number  of  adjacent  estimates 

The  larger  the  number  of  adjacent  estimates  that  are  averaged,  the  greater  the  accuracy, 
but  not  without  a  price.  The  increase  in  effective  bandwidth  also  lowers  the  frequency 
resolution. 

Another  method  of  improving  the  accuracy  of  a  spectral  estimate  is  the  application 
of  a  tapered  data  window.  This  function  modifies  the  input  data  record  by  weighting  the 
middle  of  the  digital  information  more  heavily  than  the  two  ends.  The  overall  effect  is  a 
smoothing  of  the  input  data  with  a  corresponding  increase  in  the  sharpness  of  the  spectral 
window.  The  increase  in  sharpness  comes  at  the  expense  of  a  slightly  larger  bandwidth 
and  lower  frequency  resolution.  Figure  7  shows  one  type  of  a  tapered  data  window  called 
a  Hamming  window. 


Figure  7. 

Hamming  Window 
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The  following  steps,  as  given  by  Newland  [12],  summarize  the  procedure  for  aevei- 
oping  spectra  from  digital  input  data: 

1 .  Estimate  the  frequency  range  of  interest. 

2.  Choose  an  appropriate  sampling  interval  A  to  insure  that  the  Nyquist  fre¬ 
quency  exceeds  the  frequency  of  interest.  Filter  the  input  data  to  remove  fre¬ 
quencies  higher  than  the  Nyquist  frequency. 

3.  Decide  on  the  required  accuracy  c/m. 

4.  Estimate  the  required  frequency  resolution  and  specify  the  maximum  effec¬ 
tive  bandwidth  Be. 

5.  Calculate  the  total  sampling  time  T  from  Equation  4.7. 

6.  Determine  the  number  of  data  points  needed  from  N  =  T/A.  The  total  num¬ 
ber  of  data  points  for  each  record  must  be  equal  to  a  power  of  two  to  apply 
FFT  analysis. 

7.  Determine  the  number  of  adjacent  spectral  estimates  that  must  be  averaged 
from  Equation  4.8. 

8.  Conduct  the  Fourier  analysis  using  data  taper  windows  if  desired. 

9.  Average  the  required  number  of  adjacent  spectral  estimates  to  obtain  the 
smooth  spectrum. 

4.2  SPECTRUM  DEVELOPMENT 

The  spectra  in  Appendix  B  were  generated  using  the  theory  and  procedure  outlined 
in  Section  4.1.  Table  II  lists  the  particular  parameters  used  in  the  Fourier  analysis. 
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Table  II. 

Fourier  Analysis  Parameters 


Frequency  Range 

0-500  Hz 

Nyquist  Frequency  fa 

1000  Hz 

Sampling  Frequency 

2000  Hz 

Sampling  Interval  A 

0.0005  sec 

Number  of  Sample  Points  N 

4096 

Total  Sample  Time  T 

2.0475  sec 

Number  of  Estimates  Averaged 

5 

Effective  Bandwidth  Be 

2.442  Hz 

Frequency  Resolution 

4.884  Hz 

Accuracy  c/m 

0.4472 

The  data  for  the  spectra  was  gathered  by  surveying  the  area  just  forward  of  the  foil, 
the  incoming  boundary  layer,  and  the  free  stream.  Measurements  were  made  at  tunnel 
speeds  of  10  and  20  ft/sec.  This  equates  to  Reynolds  numbers  of  8.7xl05  and  1.7xl06 
based  on  the  chord  length  of  the  foil.  If  Reynolds  numbers  are  calculated  using  the 
momentum  thickness  of  the  theoretical  boundary  layer,  as  in  previous  studies,  the  results 
are  5330  and  9280  respectively.  Attempts  were  also  made  to  collect  data  at  a  tunnel 
speed  of  25  ft/sec,  but  tracker  lock-on  could  not  be  maintained. 

The  coordinate  system  shown  in  Figures  1  and  2  was  used  as  a  reference  for  the 
data  collection.  The  origin  of  the  coordinate  system  was  at  the  junction  of  the  foil  and  the 
splitter  plate  as  indicated  in  Figure  2.  The  location  of  each  measurement  point  was  given 
in  millimeters  to  allow  use  of  the  traverse  table  software.  The  area  forward  of  the  sail 
was  surveyed  on  the  plane  of  symmetry  using  a  two  dimensional  grid  of  six  horizontal 
points  and  three  vertical  points.  The  size  of  the  grid  was  set  to  incorporate  the  significant 
portion  of  the  junction  vortex.  The  boundary  layer  was  sampled  at  a  position  255  mm 
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upstream  of  the  foil,  the  same  location  used  to  measure  the  boundary  layer  profiles  given 
in  Appendix  A.  The  free  stream  data  was  collected  in  the  test  section  with  only  the  split¬ 
ter  plate  installed. 

The  data  collected  by  the  analog-to-digital  acquisition  system  was  placed  on  floppy 
disks  for  off-line  analysis.  The  first  step  of  the  analysis  was  to  put  each  data  file  through 
a  digital  4th-order  Butterworth  bandpass  filter.  For  most  of  the  spectra  in  Appendix  B, 
the  bandpass  filter  was  set  at  10-1000  Hz.  The  upper  limit  was  set  to  match  the  setting  of 
the  analog  low-pass  filter  of  the  signal  conditioner  in  the  acquisition  system.  Once  again, 
the  value  of  this  cutoff  was  determined  from  the  Nyquist  frequency.  The  lower  limit  was 
chosen  to  remove  low  frequency  fluctuations  in  the  voltage  signal  caused  by  fluctuations 
in  the  tracker  data  rate.  Figures  8  and  9  show  the  effect  of  the  bandpass  filtering. 
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Figure  8. 

Unfiltered  Data 
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Figure  9. 

Data  after  10-1000  Hz  Bandpass  Filter 

Figure  8  was  produced  using  the  same  data  as  Figure  5,  except  the  data  has  been  time 
averaged  to  reduce  the  clutter  of  the  display.  Note  the  fluctuation  in  the  voltage  starting 
at  approximately  0.5  seconds.  This  was  caused  by  a  drop  in  the  tracker  data  rate.  As  can 
be  seen  in  Figure  9,  the  bandpass  filter  removed  the  fluctuation  providing  smooth  data  for 
the  spectral  analysis.  Also  note  that  the  mean  value  or  0  Hz  component  of  the  voltage 
was  removed.  This  was  of  no  consequence  since  the  spectral  analysis  software  would 
have  removed  this  component  prior  to  generating  the  spectrum. 

The  filtered  data  was  then  used  to  generate  the  power  spectral  density  (PSD)  for 
each  measurement  point.  In  all  cases,  the  spectra  were  normalized  to  a  maximum  value 
of  1.0.  Because  of  this  normalization,  some  of  the  data  files  were  filtered  using  a  50-250 
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Hz  bandpass  filter  and  then  analyzed  to  enhance  the  spectral  peaks  in  the  100  to  200  Hz 
range.  These  spectra  can  also  be  seen  in  Appendix  B.  A  Hamming  tapered  data  window 
was  applied  to  all  data  files  during  the  analysis  to  enhance  and  sharpen  the  spectral  peaks. 
A  typical  output  of  the  spectral  analysis  is  shown  in  Figure  10. 
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Figure  10. 

Power  Spectrum  using  Hamming  Window 

The  final  step  necessary  to  produce  the  smooth  spectra  in  Appendix  B  was  to  per¬ 
form  the  averages  of  the  adjacent  spectral  estimates.  Five  adjacent  values  were  averaged 
producing  the  accuracy  given  in  Table  II.  Figure  1 1  depicts  a  final  smooth  output 
spectrum. 
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Smooth  Power  Spectrum 
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5  RESULTS 


Appendix  A  contains  the  measured  profiles  of  the  approach  boundary  layers  for  tunnel 
speeds  of  10,  20,  and  25  ft/sec.  As  stated  before,  these  profiles  were  measured  at  a  position 
255  mm  upstream  of  the  leading  edge  of  the  foil  on  the  plane  of  symmetry.  The  purpose  of 
these  measurements  was  to  insure  that  a  turbulent  boundary  layer  was  properly  developing 
along  the  splitter  plate.  The  flow  was  tripped  just  downstream  of  the  leading  edge  of  the 
plate  to  induce  this  development  (see  Section  2.2). 

The  curves  shown  in  Figures  12,  13,  and  14  are  the  profiles  predicted  theoretically 
assuming  a  l/7th  power  law  for  the  velocity  distribution  and  the  measured  thickness  of  the 
boundary  layer.  Each  experimental  data  point  was  determined  by  averaging  the  values  of 
1024  samples  taken  over  a  200  sec  time  period.  The  data  shows  a  great  deal  of  scatter,  espe¬ 
cially  at  10  ft/sec,  much  more  than  would  be  expected.  The  higher  the  flow  rate,  the  closer 
the  data  approaches  the  predicted  profile.  In  all  cases,  the  measured  data  generally  follows 
the  predicted  velocity  profile;  therefore,  the  boundary  layers  were  accepted  as  satisfactory  for 
the  remainder  of  the  study. 

The  first  group  of  spectra  in  Appendix  B  are  those  for  the  free  stream  at  10,  20,  and  25 
ft/sec.  At  least  two  measurements  were  made  and  analyzed  at  each  speed  to  determine  the 
reproducibility  of  the  spectra  from  the  digitized  analog  data.  These  spectra  also  serve  as  a 
baseline  to  which  the  spectra  in  the  boundary  layer  and  junction  region  can  be  compared. 

The  consistency  of  the  free  stream  spectra  at  10  and  20  ft/sec  is  limited.  In  both  cases, 
the  approximate  size  and  number  of  frequency  peaks  are  similar  between  measurements,  but 
the  location  or  frequencies  of  the  peaks  appear  shifted.  The  magnitude  and  direction  of  the 
shift  varies  depending  on  the  frequency,  possibly  indicating  that  the  source  of  these  various 
frequencies  is  not  common  to  all.  The  majority  of  the  peaks  do  not  shift  by  more  than  50  Hz, 
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but  there  are  some  exceptions.  If  the  first  spectrum  from  10  ft/sec  is  neglected,  the  remaining 
two  are  almost  exactly  alike  with  the  exception  of  the  frequency  shifts,  especially  in  the  500 
Hz  to  1000  Hz  region. 

The  collection  of  spectra  generated  using  the  data  from  the  boundary  layer  and  junction 
region  are  extremely  disappointing.  They  appear  to  contain  very  little,  if  any,  quantitative 
information.  In  view  of  the  inconsistency  of  the  free  stream  spectra,  no  particular  frequen¬ 
cies  can  be  isolated  with  any  degree  of  accuracy.  All  attempts  to  correlate  a  spectrum  within 
the  junction  region  with  other  spectra  in  the  region,  the  boundary  layer,  or  the  free  stream 
prove  unsuccessful.  No  distinct  frequencies  or  frequency  pattern  can  be  identified  and 
directly  attributed  to  the  flow  within  the  region. 

A  small  amount  of  qualitative  information  may  be  extracted  from  the  spectra  despite 
their  limitations.  Near  the  plate  (Z  =  1mm),  at  a  tunnel  speed  of  10  ft/sec,  the  spectrum  for 
X  =  5mm  indicates  a  very  insignificant  distribution  of  energy  in  the  frequency  range  of  inter¬ 
est  (0-500  Hz).  This  spectral  energy  distribution  changes  significantly,  however,  as  the  posi¬ 
tion  of  measurement  point  moves  away  from  the  leading  edge  of  the  foil  along  the  j  n<  of 
symmetry.  The  spectrum  from  data  taken  at  X  =  10mm  shows  a  large  increase  in  all  spectral 
components  producing  numerous  frequency  spikes.  Based  on  the  predictions  of  previous  stu¬ 
dies  and  knowing  the  location  of  the  flow’s  stagnation  point  from  mean  velocity  measure¬ 
ments,  this  position  is  close  to  the  core  of  the  junction  vortex.  As  the  measurement  point 
continues  to  move  away  from  the  leading  edge,  the  trend  appears  to  reverse  resulting  in  lower 
average  spectral  levels  as  the  spectra  approach  the  levels  found  in  the  boundary  layer. 

It  is  difficult  to  confirm  the  above  trend  in  spectral  energy  levels  from  the  data  gathered 
at  20  ft/sec.  In  general,  the  levels  measured  in  the  junction  region  at  this  tunnel  speed  are 
significantly  lower  than  those  recorded  at  10  ft/sec. 


Page  34 


The  velocity  probability  distributions  in  Appendix  C  are  much  more  conclusive  in  pro¬ 
viding  useful  information  on  the  nature  of  the  junction  flow.  These  distributions  were 
created  using  the  new  digital  flow  analyzers  in  conjunction  with  the  laser  doppler  system. 
The  low  data  rate,  which  made  these  analyzers  impractical  for  the  spectra  development,  was 
not  a  factor  in  this  portion  of  the  study.  The  IFA’s  inherent  accuracy,  however,  was  crucial 
for  this  part  of  the  investigation.  The  distributions  in  Appendix  C  clearly  show  the  bi-modal 
characteristic  in  the  structure  of  the  junction  root  vortex. 

The  velocity  probability  distribution  of  the  free  stream  at  10  ft/sec  is  exactly  what 
would  be  expected.  The  shape  is  Gaussian  with  a  very  small  standard  deviation.  The  slight 
unsymmetrical  nature  of  the  curve  about  the  mean  is  caused  by  the  input  filter  settings  of  the 
IFA.  The  turbulence  intensity  calculated  from  this  distribution  is  less  than  0.5  percent  when 
adjusted  for  the  phase  noise  of  the  laser  signal.  The  distribution  in  the  boundary  layer  is  also 
a  Gaussian  curve  with  a  larger  standard  deviation  caused  by  the  increased  turbulence  in  this 
region. 

The  velocity  probability  distributions  within  the  junction  region  are  unlike  those  of  the 
free  stream  and  boundary  layer.  In  all  cases,  the  distributions  appear  to  be  the  sum  of  two 
Gaussian  distributions  displaying  two  distinct  peaks.  The  relative  strength  of  the  two  peaks, 
as  expected,  is  dependent  on  the  distance  away  from  the  leading  edge  of  the  foil.  By  follow¬ 
ing  the  progression  of  figures  in  Appendix  C,  it  is  obvious  that  the  larger  peak  shifts  from  a 
negative  velocity  behind  the  stagnation  point  to  a  positive  velocity  as  the  stagnation  point  is 
crossed.  The  presence  of  two  peaks  is  indicative  of  two  distinct  flow  velocities  existing  over 
time  at  the  measurement  point  confirming  a  bi-modal  structure. 
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The  velocity  probability  distribution  results  obtained  at  a  tunnel  speed  of  20  ft/sec  are 
the  same  as  those  acquired  at  10  ft/sec.  The  bi-modal  characteristic  of  the  flow  structure  is 
even  more  pronounced  at  this  tunnel  speed.  The  two  distinct  peaks  are  clearly  visible  at  X  = 
5mm  and  X  =  10mm. 
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6  CONCLUSIONS  AND  RECOMMENDATIONS 


The  conclusions  that  can  be  drawn  from  this  experimental  study  of  the  junction  root 
vortex  are  not  as  extensive  as  originally  expected.  The  results  of  the  velocity  probability  dis¬ 
tributions  clearly  confirm  the  presence  of  the  bi-modal  characteristic  within  the  structure  of 
the  junction  vortex  as  stated  by  Simpson.  The  spectral  studies,  however,  provide  no  clear 
information  about  the  flow  except  to  indicate  an  increase  in  the  spectral  energy  distribution  in 
the  vicinity  of  the  vortex.  The  frequency  of  fluctuation  of  the  bi-modal  structure  or  any  other 
distinctive  frequency  pattern  caused  by  the  flow  cannot  be  identified  from  the  spectra 
obtained  in  this  investigation. 

The  ability  to  use  data  from  MIT’s  laser  doppler  anemometry  system  to  obtain  spectral 
information  is  still  questionable  and  requires  further  study.  The  accuracy  of  the  data  gath¬ 
ered  using  the  IFA’s  could  significantly  improve  the  quality  of  the  spectral  output  if  suffi¬ 
cient  data  rates  could  be  obtained. 

The  following  is  a  list  of  recommended  procedures  for  further  investigation.  Some  of 
these  procedures  relate  directly  to  the  investigation  of  the  junction  vortex,  while  others 
broadly  apply  to  the  use  of  the  laser  doppler  system  for  time-frequency  studies: 

1.  Conduct  on-line  frequency  analysis  of  the  laser  doppler  signal  using  a  spec¬ 
trum  analyzer  and  the  analog  output  of  a  signal  processor9.  The  sample-and-hold 
nature  of  both  the  tracker  and  IFA  must  be  taken  into  account.  A  sufficient  data 
rate,  depending  on  frequencies  of  interest,  must  be  achievable  to  assume  a  contin¬ 
uous  input  signal  to  the  analyzer. 

2.  Improve  the  data  rates  obtainable  with  the  IFA  to  allow  the  collection  of 


9  Although  the  primary  output  of  an  IFA  is  digital  information,  an  analog  output  is  available  from  a  connection 
on  the  back  panel  of  the  processor. 
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accurate  digital  information  for  off-line  spectral  analysis  within  a  useful  fre¬ 
quency  range. 

3.  If  sufficient  data  rates  can  be  obtained  from  the  IFA’s  on  both  the  horizontal 
channel  and  fiber  optic  (third  component)  channel,  attempt  to  cross-correlate 
fluctuations  in  the  junction  region  with  fluctuations  in  the  free  stream  or  bound¬ 
ary  layer. 

4.  Improve  the  accuracy  of  the  spectra  by  averaging  a  larger  number  of  adjacent 
spectral  estimates.  A  a/m  of  ten  percent  can  be  obtained  by  averaging  100  adja¬ 
cent  /alues. 

5.  Average  a  large  number  of  spectra  generated  from  sets  of  data  taken  at  the 
same  point  over  time.  This  may  help  to  reduce  the  fluctuations  in  the  spectra,  as 
indicated  in  the  free  stream  cases,  allowing  only  significant  frequency  compo¬ 
nents  to  emerge. 
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APPENDIX  A 


BOUNDARY  LAYER  PROFILES 
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Figure  12. 

Boundary  Layer  Profile  (10  ftysec) 
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Boundary  Layer  Profile  (20  ftysec) 
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APPENDIX  B. 


POWER  SPECTRAL  DENSITIES 
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Figure  16. 

Free  Stream,  10  ft/ sec 
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Figure  19. 

Free  Stream,  10  ft/sec 
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Figure  22. 

Free  Stream,  20  ft/sec 
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Figure  24. 

Free  Stream,  20  ft/sec 


0.9 


Page  49 


Pr®qu«ncy  (KHz) 


Figure  25. 

X=5mm,  Z=lmm,  10  ft/sec 
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Figure  26. 

X=5mm,  Z=lmm,  10  ft/sec 
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Figure  29. 

X=10mm,  Z=lmm,  10  ft/sec 
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Figure  31. 

X=15mm,  Z=lmm,  10  ft/sec 
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Figure  32. 

X=15mm,  Z=lmm,  10  ft/sec 
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Figure  33. 

X=15mm,  Z=lmm,  10  ft/sec 
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Figure  35. 

X=25mm,  Z=lmm,  10  ft/sec 
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Figure  36. 

X=25mm,  Z=lmm,  10  ft/sec 
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Figure  38. 

X=5mm,  Z=8mm,  10  ft/sec 
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Figure  41. 

X=15mm,  Z=8mm,  10  ft/sec 
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Figure  42. 

X=15mm,  Z=8mm,  10  ft/sec 
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Figure  43. 

X=25mm,  Z=8mm,  10  ft/sec 


Figure  44. 

X=25mm,  Z=8mm,  10  ft/sec 
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Figure  46. 

X=5mm,  Z=2.8mm,  10  ft/sec 
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Figure  50. 

X=15mm,  Z=2.8mm,  10  ft/sec 
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Figure  5 1 . 

X=15mm,  Z=2.8mm,  10  ft/sec 
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Figure  52. 

X=25mm,  Z=2.8mm,  10  ft/sec 
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Figure  53. 

X=25mm,  Z-2.8mm,  10  ft/ sec 
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Figure  56. 

Boundary  Layer,  10  ft/sec 
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Figure  57. 

Boundary  Layer,  10  ft/ sec 
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Figure  59. 

X=5mm,  Z=lmm.  20  ft/sec 
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Figure  61. 

X=10mm,  Z=lmm,  20  ft/sec 
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Figure  62. 

X=15mm,  Z=lmm,  20  ft/ sec 
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Figure  63. 

X=lSmm,  Z=lmm,  20  ft/«ec 
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Figure  65. 

X=25mm,  Z=1  mm,  20  fl  sec 
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Figure  67. 

X=5mm,  Z=7mm,  20  ft/sec 
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Figure  69. 

X=15mm,  Z=7mm,  20  ft/sec 


Page  77 


PSD 


0.5  0.7  0  9 

Fr®qu«ncy  (KHz) 


Figure  71. 

X=25mm,  Z=7mm,  20  ft/ sec 
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Figure  73. 

X=5mm,  Z=2.1mm,  20  ft/sec 
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Figure  74. 

X=15mm,  Z=2.1mm,  20  ft/sec 
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Figure  75. 

X=15mm,  Z=2.1mm,  20  ft/ sec 
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Figure  76. 

X=25mm,  Z=2.1mm,  20  ft/sec 
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Figure  77. 

X=25mm,  Z=2.1mm,  20  ft/ sec 
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Figure  78. 

Boundary  Layer,  20  ft/sec 
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Figure  79. 

Boundary  Layer,  20  ft/sec 
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Figure  80. 

Free  Stream,  10  ft/sec 

The  PSD  ic  Gaussian  with  a  small  standard  deviation. 

The  non-symmetr>  o  due  to  input  filter  settings  of  the  signal  processor. 


Page  8 


VELOCITY  (M/S) 


Figure  85. 

X=15tnm,  Z=lmm,  10  ft/  sec 

The  dominant  velocity  shifts  to  positive  flow. 
Bimodal  structure  is  still  evident. 
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Figure  89. 

X = 1 .5  m  m,  Z=lmm,  20  ft/sec 

Double  peak  is  already  starting  to  become  evident  very  close  to  foil. 
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Figure  93 

X=20mm,  Z=lmm,  20  ft/sec 
approaches  the  shape  of  the  boundary  layer  distribution 
Tne  double  peak  is  still  slightly  evident. 


